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We present an extinction map of the Large Magellanic Cloud (LMC), using 
204,502 stars from the Two Micron All Sky Survey point source catalog. We first 
use the NICE method to determine the reddening distribution, E(H — K) and 
E(J — H), which we compare to the HI distribution to find a near-infrared reddening 
law of E(J - H)/E(H - K) = 1.20 ± 0.04. A visual extinction map (~ 6° x 6°) 
of the LMC is created using the NICER method; at 4 arcmin resolution, a mean 
value of Ay = 0.38 mag is found. We derive the LMC CO-to-H2 conversion fac- 
tor, Xlmc, independent of assumptions about the virialization of giant molecular 
clouds, by comparing the NICER extinction map with NANTEN 12 CO observa- 
tions. In regions where Ay > 1 mag and 12 CO emission is > 2 K km s _1 , we 
measure Almc — 9.3 ±0.4 x 10 20 cm -2 (K km s -1 ) -1 . In the same regions, the LMC 
contains a total molecular mass of (4.5 ± 0.2) x 10 7 M . 

Subject headings: dust, extinction — galaxies: ISM — ISM: abundances — ISM: clouds 
— Magellanic Clouds — techniques: photometric 



1. Introduction 

Studying the properties of giant molecular 
clouds (GMCs), the sites of most star for- 
mation, in a variety of extragalactic environ- 
ments provides insight into stellar evolution as 
a whole. The Large Magellanic Cloud (LMC) 
is a good laboratory for this type of study be- 
cause of its proximity (D ~ 50 kpc) and its 
nearly face-on orientation in the sky. 

Determining the amount and distribution 
of molecular hydrogen (H 2 ), the dominant 
molecule (~ 99.99% by number) in GMCs, 
is useful for the study of the initial conditions 
of star formation. Yet due to the lack of a 
permanent electric dipole moment, H2 is un- 
detectable at the low temperatures (~ 10 K) 
of molecular clouds, and so a variety of other 
tracers have been used to identify them and 
to study their density distribution. Two com- 
mon methods are star counts and color excess, 
which trace molecular gas by measuring the 
dust distribution of GMCs. The spatial lo- 



cation of gas and dust are highly correlated, 
as the well-known gas-to-dust relationship at- 
tests, because H 2 forms most efficiently in 
dense gas where dust grains provide surfaces 
onto which H I can collide and be converted 
into molecular form. The star counts method 
measures extinction by comparing the num- 
ber of stars seen behind a dark cloud along 
a given light of sight with the number in a 
control field off the clouds and assumed to be 
dust-free. The color excess method measures 
reddening by exploiting the decrease in the 
extinction with increasing wavelength. 

Since the 1960s, a number of authors have 
used the color excess method to map the red- 
dening and inferred dust distribution of the 
LMC, (e.g., Feast et al. 1960; Isserstedt 1975; 
Grieve & Madore 1986). But a combination 
of questionable assumptions — for instance, as- 
suming Galactic intrinsic colors, small data 
sets, limited spatial coverage, and not account- 
ing for the Galactic foreground contribution to 
reddening — resulted in significant uncertain- 
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ties. 

Advances in recent years have made it pos- 
sible to possible to bypass most of these prob- 
lems. Namely, the advent of infrared (IR) ar- 
ray detectors, such as those used in the Two 
Micron All-Sky Survey (2MASS), have made 
possible the simultaneous observations of hun- 
dreds of sources in multiple wavelength bands. 
These detectors have the further advantage 
that IR light suffers much less extinction than 
the visible and ultraviolet, so molecular clouds 
can be probed more deeply. In the Galaxy, 
for instance, Ak ~ 0.1 Ay. The Near Infrared 
Color Exces s (NIC E) technique developed by 
Lada et al. (119941 ) exploits the advantages of- 
fered by these large-scale arrays and has proved 
to be a powerful way of mapping extinction 
through molecular clouds. In this method, in- 
frared color excess measurements are made by 
inferring the intrinsic color, typically (H — K), 
of target stars from an unreddened control 
field. The NICER (NICE Revised) method of 
Lombardi & Alves (2001) generalizes the for- 
mer by using information from all available in- 
dependent colors — for instance, (J — H) and 
(H — K) — to map extinction. 

Radio observations of interstellar carbon 
monoxide (CO) h ave a lso been used to trace 
H 2 . Cohen et al. (119881 ) were the first to fully 
survey the LMC in the J = 1 —>■ rotational 
transition of CO. Yet, at their limited resolu- 
tion (8/8), they were able to resolve only the 
largest cloud complex es (~ 140 pc). More re- 
cently, Mizuno et al. (l200lh completed a com- 
prehensive high-resolution (~ 40 pc) CO sur- 
vey using the 4 m NANTEN telescope. Both 
groups of authors made calculations of the CO- 
to-H2 conversion factor, or the "A-factor," of 
the LMC. 

An outstanding question concerning GMCs 
is how metallicity affects the A-factor. This 
quantity is of particular interest in the LMC 
since the galaxy has a lower metallicty (~ 1/4; 
Dufour 1984) and a higher gas-to-dust ratio 
(~ 1.7; Gordon et al. 2003) than the MW. It 
is conjectured that the relatively less-abundant 
CO in the LMC is more supceptible to photo- 



dissociation by UV radiation, which can pen- 
etrate more deeply into the less-dusty clouds 
(e.g., McKee 1989). This argues for X LMC > 
Amw- 

Yet there number of problems with 

determining the A-factor soley from CO ob- 
servations. Potentially unrealistic assump- 
tions regarding the physical properties of the 
molecular clouds — for instance, that they are 
virialized — must be made, since there is not an 
independent measure of A(H 2 ). 

It is the main goal of this paper to provide 
such an independent measure of the molecu- 
lar content in the LMC, one that makes no 
assumption about the structure, dynamics, or 
virialization of GMCs. We do this by first cre- 
ating an extinction map from 2MASS data, in 
one of the first applications of the NICE(R) 
techniques to an extragalactic system. A num- 
ber of other useful results are obtained in the 
process, including the NIR extinction law, the 
extinction distribution, and the A-factor of the 
LMC. In § 2 we describe data selection and the 
other considerations that go into determining 
extragalactic reddening. The results, including 
the extinction map and A-factor, are presented 
in § 3. Conclusions are presented in § 4. 

2. Data Selection 

2.1. 2MASS Data 

Our reddening maps of the LMC are derived 
using JHK near-IR photometry obtained from 
the 2MASS Point Source Catalog^]. For uncon- 
fused sources with Galactic latitude |6| > 10°, 
the limiting magnitudes of the 2MASS catalog 
at the 10-cr level are 15.8, 15.1, and 14.3 mag in 
J (1.24 /mi), H (1.66 /on), and K s (2.16 /mi), 
respectively. Possible contaminants, such as 
image artifacts caused by bright stars, con- 
fused sources, solar system objects, and other 
artifacts are eliminatecH. Only data having 



1 http://irsa. ipac.caltech.edu/applications/Gator/ 
2 Explanatory Supplement to the 

2MA SS All- Sky Data Release at 
|http:/ / www.ipac.caltecli.edu /2mass/releases/ all- 
sky/ doc/explsup.html 
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photometric uncertainty < 0.10 mag, corre- 
sponding to S/N > 10, are selected. Our data 
cover the central 6° x 6° region of the LMC, 
centered on (a, 5) = (5 h 20 m , -69°). This 
choice of range was b ased o n an examination 
of the Mizuno et al. (2001) CO map and the 
Staveley-Smith et al. (120031 ) H I map, which 
reveals the extent of molecular and H I emis- 
sion, respectively. This initial sample contains 
329,702 sources. 

The 2MASS resolution, 2", corresponds to 
a physical size of ~ 0.5 pc at 50 kpc, and so 
individual LMC stars should be well-resolved. 



Nikolaev & Weinberg (120001 ). treating a simi- 
lar set of LMC data drawn from 2MASS, show 
that confusion is not a problem at this resolu- 
tion. 

Our preliminary data set contains not only 
LMC field stars but Galactic foreground point 
sources. In order to get a quantitative esti- 
mate of the MW foreground stellar number 
density, we obtain 2MASS data for a reference 
field in the same range of Galactic latitude as 
the LMC field but removed in Galactic longi- 
tude by —30°. The angular number densities 
of the LMC field (prior to data reduction) and 
MW field are 4.65 x 10 3 deg" 2 and 1.68 x 10 3 
deg -2 , respectively. Not accounting for the 
foreground would dilute the measured extinc- 
tion (Ay) and increase the uncertainty, espe- 
cially in regions of high Ay. In what follows 
we use the color-color diagra m, in conjunction 
with Nikolaev & Weinberg's (120001 ) analysis of 
the LMC color-magnitude diagram, to correct 
for this. 

2.2. Data Reduction 

To derive the LMC gas-to-dust ratio (see 
§ 13.11) and X-f actor (§ 13.31) . we generate color 
excess and visual extinction maps for compari- 
son with the galaxy's N(H I) and CO distri- 
butions, respectively. We use Lada et al.'s 
(1994) NICE method to generate the color ex- 
cess maps, E(J — H) and E(H — K). For the 
extinction maps, we use the NICER method 
of Lombardi & Alves (2001), which general- 
izes Lada's formulation for multi-band obser- 



vations. In these techniques, the intrinsic color 
of all stars observed toward a target clouds 
is taken to be the mean color of a nearby, 
unreddened control field. An extinction map 
is generated by averaging the irregularly sam- 
pled data into pixels, across each of which the 
extinction is assumed to be uniform. Appli- 
cation of these techniques, originally applied 
to Galactic clouds, assumes that all stars ob- 
served toward the cloud (i.e., the field stars) 
are background stars; and the field and con- 
trol stars are homogeneous (i.e., drawn from 
the same parent population). Extending these 
techniques to extragalactic sources raises spe- 
cial concerns as to how well these conditions 
can be fulfilled. 

Before creating the extinction map, a num- 
ber of reductions and corrections are per- 
formed on the raw data, with particular care 
taken to address those assumptions just men- 
tioned. 

2.2.1. Completeness 

The observed completeness limits toward 
the LMC are brighter than the nominal com- 
pleteness limits stated in the 2MASS catalog 
(Figured]). Towards regions of higher extinc- 
tion, intrinsically fainter stars will escape de- 
tection before intrinsically brighter stars, so 
the overall extinction would be underestimated 
if all stars were included in the measurements 
(Oestreicher & Schmidt-Kaler 1996). Because 
only the minimally extincted faint stars would 
be able to be measured, including them in our 
sample would drive down the overall mean ex- 
tinction. However, in dealing with apparent, 
not absolute, magnitudes, the relative faint- 
ness of a given star may be due to extinc- 
tion or its intrinsic brightness. But statisti- 
cally, we do not expect a particular popula- 
tion (intrinsically faint or bright) to be more 
likely to lie behind molecular clouds. With 
our large data set, then, it is reasonable to as- 
sume that most observed faint stars (i.e., those 
below the observed completeness limit) would 
bias measurements to lower values of extinc- 
tion if included. Thus, based on our complete- 
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Fig. 1. — Luminosity functions of our sam- 
ple. The observed completeness limits for J, 
H, and K, were estimated 15.0, 14.4, and 14.3 
mag, respectively. The 2MASS completeness 
limits (15.8, 15.1, and 14.3 mag) are indicated 
with vertical lines. 

ness estimates, we used only stars with bright- 
ness greater than J < 15.0, H < 14.4, and 
K < 14.3 mag for subsequent analaysis. This 
and the following data selections are summa- 
rized in Table [TJ 

We tested these assumptions in several ways. 
First, we confirmed that the use of a constant 
magnitude limit in each band is appropriate, 
i.e., that it does not depend on column density. 
We calculated the limits for different redden- 
ing intervals and found that the values stated 
above remained the same. Second, we mea- 
sured the extinction for those stars falling out- 
side our stated limits (i.e., those fainter than 
J > 15.0, H > 14.4, and K > 14.3 mag). 
Their measured mean reddening of -0.44 mag 
is unphysical and their distribution is consis- 
tant with an error population. Indeed, when 
we went ahead and created a map including 
these stars, (i.e., without having imposed a 
completeness limit), the global mean extinc- 
tion was lowered significantly (by nearly a fac- 



tor of 2). Finally, after having created the ex- 
tinction map only using stars brighter than the 
completeness estimates, we found that the re- 
maining fainter stars have overall lower redden- 
ings consistant with an error population (i.e., a 
noise population with a Gaussian distribution 
centered around zero). 

2.2.2. Foreground and non-interstellar red- 
dening 

The Galactic foreground contaminates the 
2MASS LMC field in two ways: with H I gas 
and point sources. Within the LMC itself there 
is foreground contamination, as well as red- 
dening of non-interstellar origin, including that 
caused by protostars and obscured aymptotic 
giant branch (AGB) stars. 

Galactic H I foreground extinction is not 
negligible but previous authors, usually for 
lack of data, adopted uniform foreground red- 
dening in their derivation of LMC extinction 
maps. In the Parkes multibeam H I survey of 
the LMC, Staveley-Smith et al. (2003) found a 
mean Galactic reddening of (E(B — V)) = 0.06 
mag over the LMC disk, varying from 0.01 
to 0.14 mag. Staveley-Smith generously pro- 
vided us with their H I data, which we used 
to correct for Galactic reddening for stars in 
the LMC field. For every 2 MASS source, 
foreground extinction in each band was cal- 
culated at that position using the MW ex- 
tinction law of Schlegel, Finkbeiner & Davis 
fll998h (Aj/Ay = 0.276, A H /A Y = 0.176, 
A K /A V = 0.112), an d the MW gas-to-dust ra- 
tio of Bohlin et al. fll978f ) {(3 = iV(H 1)/A y = 
1.87 x 10 21 mag -1 cm -2 ). These values were 
then subtracted from the observed 2MASS 
magnitudes: 
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While Galactic H I is easily identified given 
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its distinct velocity compared to LMC gas, the 
correction for Galactic point source contami- 
nation is more complicated. Furthermore, we 
wish to eliminate from our sample objects in 
the LMC whose reddening is circumstellar in 
origin. The color-magnitude diagram (CMD) 
is a useful tool to correct for both effects. Anal- 
ysis of the CMD allows us to eliminate sources, 
on a probablistic basis, that are unlikely to 
have undergone reddening from LMC interstel- 
lar dust. 



a -69 




R.A. (J2000.0) 

Fig. 2. — Spatial distribution of the selected 
sources, in pixels of 4' x 4', with average sur- 
face density of Elmc = 4.76 x 10 3 stars deg -2 . 
Counterclockwise starting from the northeast, 
the outlined control regions have an average of 
£ cont . = 1.48, 1.57, 1.98, and 1.02 x 10 3 stars 
deg -2 . The color levels cyan, green, yellow, or- 
ange, and red correspond to 0.90, 1.80, 3.60, 
7.20, and 14.40 xlO 3 stars deg -2 , respectively. 



One can use the CMD (K vs. J-K) to 
identify distinct stellar populations (Figure [3]). 
This is the main focus of Nikolaev & Wein- 
berg fl2000h who, by matching features of the 
2MASS CMD to colors of known populations 
from the literature, discuss which populations 
are the greatest contributors to the different 
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Fig. 3. — Color- magnitude diagram of the 
LMC field. The contours are logarithmically 
spaced, from 2 to 3.4, by 0.4. We selected data 
for the extinction maps from the enclosed re- 
gion, which corresponds to regions E, F, G, H, 
J, L and part of D of Nikolaev & Weinberg 
(2000). 

regions. We use their regions E, F, G, H, J, 
K, L, and part of D, to create the reddening 
maps. These 3X6 clS correspond to LMC popu- 
lations, primarily of the giant branch, with in- 
significant Galactic foreground contamination. 
We eliminate regions A, B, C, and I, which 
contain roughly 15, 80, 80, and 55 percent 
foreground, respectively. Region K, which is 
free of foreground contamination, but consists 
of dusty asymptotic giant branch stars whose 
large J — K colors are due to circumstellar 
reddening, is also eliminated. We confirmed 
that these highly reddened stars are not pref- 
erentially located in regions of high molecu- 
lar gas by plotting their spatial locations. We 
also performed the reverse operation. That 
i s, ba sed on the CO map of Mizuno et al. 



( 120011 ). we located points on the color-color di- 
agram and CMD that fall within the contours 
of known molecular clouds. In all, there are 
only 74 stars in the extended "finger" of the 
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Fig. 4. — Color-color diagram of the LMC 
field before data reduction. The contours are 
logarithmically spaced, from 2 to 3.8, by 0.4. 
Overplotted are the color sequences for Galac- 
tic dwarfs and giants from Koornneef (1983). 
The reddening vector, based on relations from 
Koornneef (1982), is drawn for E(B-V) = 1.0 
mag. 

color-color diagram corresponding to m olecu- 
lar clouds in the Mizuno et al. (120011 ) map 
(see Figure H]). After these selections, we are 
left with 204,502 sources, corresponding to an 
average number density of 4.76 x 10 3 deg -2 . 
Figures H] and [5] show the color-color diagrams 
of the LMC field before and after the data re- 
ductions presented in this section. 

A difficult problem is correcting for fore- 
ground contamination within the LMC. The 
color excess method assumes that stars ob- 
served in the line-of-sight to a cloud are back- 
ground stars. Extinction measurements will 
be diluted wherever foreground stars are in- 
cluded in the calculations. In addition, as dis- 
cussed and modeled by Lada et al. (1994), the 
decrease in measured extinctions due to fore- 
ground stars will have a particularly large im- 
pact in regions of high extinction, resulting in 
an increase in dispersion with Ay. We tried 



Fig. 5. — Color-color diagram of the 204,502 
stars in the LMC field after data reduction. 
The contours and reddening vector are as in 
Figure BJ 

to minimize the effect of foreground in the ex- 
tinction map by sigma-clipping — that is, by re- 
moving sources in a given pixel with anomolous 
values of Ay. This was one of the motivations 
behind the creation of the NICER technique 
(Lombardi & Alves 2001). Described compre- 
hensively by these authors, the sigma-clipping 
solution is justifiable when the scatter of intrin- 
sic stellar colors is small or when the column 
density of the GMCs in question is sufficiently 
high. As we show in the following subsection, 
the first condition is not met in our case. Nei- 
ther do we have any a priori reason to sigma- 
clip based on the latter condition. 

In any case, to explore the effects of this 
method on our data set, we did generate a 
number of sigma-clipped maps. Ultimately, we 
found no advantages in doing so and used a 
simple weighted mean for our final maps. 

2.2.3. Homogeneity of control and field stars 

Success with the NICE(R) technique relies 
on the mean color of the control group being 
characteristic of the intrinsic color of the field 
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Table 1: Data selection summary. 

Criteria Number of 

remaining stars 
2 MASS selection 329,702 
Completeness cutoff 3, 258,621 
Eliminate MW foreground 15 204,502 



0.12T 




H-K (mag) 



0.08 



"Completeness limits (see Figure [T]): J < 15.0 , H < 14.4, o.06 
K < 14.3 mag. | 
b Based on Nikolaev & Weinberg ( 2000h regions. 



stars. Recall, the reddening of a given star 
is the difference between the observed color 
and the intrinsic color of a particular spectral 
type. Thus, extending the NICE(R) method 
to an entire extragalactic field raises another 
concern as to whether the fraction of stars of 
a particular spectral type varies significantly 
across the field. In other words, we are con- 
cerned with the homogeneity of the field, espe- 
cially the control group. Based on the selection 
criteria in the previous section, we are confi- 
dent that our field consists primarily of giant 
branch stars, though a variety of stages in evo- 
lution are represented (e.g., helium and hydro- 
gen shell burning phases, oxygen- and carbon- 
rich phases; Nikolaev & Weinberg 2000). As 
for the control group, we would like the disper- 
sion around the mean control color to be signif- 
icantly less than the observed range of field star 
colors, since the dispersion determines the 1-a 
confidence level for extinction measurements 
(Alves et al. 1998). 

We examined the H I and CO maps of the 
LMC for minimally reddened regions and in- 
spected control groups in each of the four "cor- 
ners" of the galaxy image (Figure [2]). Not 
surprisingly, the regions do not have nar- 
row color distributions. The control groups 
are well-behaved in that they are statistically 
equivalent — the (H — K) colors to two signif- 
icant figures and the (J — H) colors to one 
significant figure. Yet, as Figure El shows, the 
width of the control distributions relative to 
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Fig. 6. — The reddening distrutions of the 
LMC field (heavy black line) and control (col- 
ored lines) groups. The orange, red, green, and 
blue lines correspond to data from the north- 
east, southeast, southwest, and northwest cor- 
ners of the LMC, as shown in Figure [2j 

that of the field stars forecasts low signal-to- 
noise in extinction measurements. The nature 
of the uncertainties is explored in detail in 
§ 13.21 where we derive the extinction maps. 
We are somewhat consoled, though, since the 
large number of stars in our data set will act to 
minimize this problem and because our main 
interests here are in global properties of the 
galaxy. We proceed by taking the average of 
the mean control group colors as the intrinsic 
colors of the field: 

[H — K) intrinsic = \\H — K)) control 

= 0.16 ±0.09 mag, (4) 

\J Hj intrinsic = \\J -^0) control 

= 0.74 ±0.11 mag. (5) 

By comparison, the mean colors of the entire 
data set are (H-K) = 0.175 ± 0.0946 and 
(J - H) = 0.772 ± 0.110 mag (see Figure E>. 
Again, this and the previous selection criteria 
are summarized in Table [TJ 
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3. Results 

3.1. NIR Extinction Law 
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Fig. 7. — LMC reddening as a function of 
hydrogen column, measured in regions where 
Icq < 1 K km s -1 . Error bars represent 
the deviation around the mean of color ex- 
cesses inside the bin (AJV(Hl)/[2. x 3.25 x 
0.05]). 
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The NICE(R) methods require knowledge of 
the NIR extinction law and assume that this 
relationship is linear. The extinction law may 
be described by extinction coefficients, fcj = 
Ay/E(m\ 1 -m^ 2 ), which relate the extinction 
to the color excess. Knowing ki enables us to 
convert color excess to Ay, which is propor- 
tional to iV(Htotai) and, hence, to iV(n 2 ). To 
determine the scaling we must assume a value 
of the LMC gas-to-dust ratio. Thus, we plot 
E{m\ 1 —m\ 2 ) versus iV(H I), and compare this 
relationship with the LMC gas-to-dust ratio es- 
timated by Gordon et al. (2003) (Figure [7]). 
Their value of iV(H I) /Ay = (3.25±0.28) x 10 21 
mag -1 cm -2 is a weighted average of all 24 
archival gas-to-dust ratio measurements. 

E(H — K) and E(J — H) for individual stars 
were determined by subtracting the intrinsic 
colors (Equations H] and respectively) from 
the observed 2MASS colors. Color excess maps 
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Fig. 8. — The reddening law measured in the 
same regions as Figure The (H — K) color 
excess is plotted clS 8b function of (J — H) 
color excess. Error bars represent the devia- 
tion around the mean of E(H — K) inside the 
bin AE(J-H) = 0.01 mag. The data are plot- 
ted function of E( J — H), as opposed to 
E(H — K) as is often the case, because in this 
circumstance E(J — H) has a higher overall 
signal-to- noise. 

were generated by binning the data into 4' 
pixel grids with 2' spacing. The color excess 
at a given location was then given by the av- 
erage, weighted by the uncertainties, of stars 
within that pixel. Next, these maps were com- 
pared to the iV(H I) map on a pixel-by-pixel 
basis. 

In Figure the data are averaged in bins of 
equal iV(H I) and fitted for the linear relation 



E(m Xl - m X2 



rriiA^ 



(HI) +6, (6) 



where m, = k~ l and Ay(R I) = iV(H I)/(2 x 
3.25 x 10 21 mag -1 cm -2 ) is the extinction due 
to H I. The H I column is divided by 2 since, on 
average, extinction along a given line-of-sight 
is caused by half the dust layer. In other words, 
only half the total column, on average, con- 
tributes to the measured reddening. The data 
are plotted in bins of Ay4 v (H I) = 0.05 mag, 
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Fig. 9. — Sigma-clipped NICER extinction 
map of the LMC at a resolution of 4 arcmin. 

and weighted least-squared fits were done over 
the range Ay(H I) = [0.1,1] mag. So as to 
get the best possible correlation between color 
excess and atomic hydrogen we masked out re- 
gions with significant molecular gas, defined as 
areas with CO emission exceeding Iqo — 1 
K km s _1 . This value is based on an exami- 



nation of the Mizuno et al. (120011 ) CO map, 
in which they define GMC complexes as reions 
exceeding 3 K km s _1 . The extinction coeffi- 
cients ki, which came from the inverse of the 
best fit slopes, were found to be 20.83 ± 0.52 
and 17.30 ±0.46 for E(H - K) and E(J-H), 
respectively. Dividing these results gives a ra- 
tio of E(J - H)/E{H - K) = 1.20 ± 0.04, 
which is consistent with the Koornneef (1982) 
reddening law, which is E(J — H) / E(H — K) = 
1.06 when converted to the 2MASS photomet- 
ric system (Carpenter 2001). 

In principle, one should be able to deter- 
mine the NIR reddening law by comparing 
E(H - K) and E(J - H) directly. And if the 
color excesses are linearly dependent, the ra- 
tio of extinctions in two different bands is con- 
stant, a key assumption of the NICE(R) tech- 
niques. We found that this direct comparison 
was not a robust way of measuring the red- 



dening law, in part, because of the lack of a 
linear relation for low values of (H — K) and 
(J — H) — a result of our data selection (see 
Figure \5§ . When we examined the reddening 
law in the same regions as the gas-to-dust ra- 
tio analysis (i.e., where Ico < 1 K km s _1 ), 
we found that the data could be fit by two 
different slopes. Figure [8] shows that the ma- 
jority of stars with little to no reddening are 
fit by a shallower slope than more highly red- 
dened stars. Below E( J — H) pa 0.05 mag, the 
data are consistent with a noise population and 
there is no robust linear fit. Above E(J—H) pa 
0.1 mag, however, the data are well approxi- 
mated by the linear relationship E(H — K) = 
(0.84 ± 0.06) E( J -H)- (0.037 ± 0.007). The 
reciprocal of the slope, 1.19, is in good agree- 
ment with the E(J - H)/E(H - K) ratio of 
1.20 we derive via the gas-to-dust ratio mea- 
surements and with the Koornneef (1982) red- 
dening law. We also repeated this analysis for 
the entire reddening maps, i.e., not restricted 
to regions of low molecular gas. The same re- 
sult was observed: more highly reddened stars 
were fit by a linear relationship in agreement 
with a reddening law with a slope of ~ 1.20. 

Finally, we note that the lack of linearity in 
the reddening law over the entire range of val- 
ues is not due to poor photometry. When we 
repeated the above analysis using stars with in- 
creasingly higher S/N ratios, we achieved sim- 
ilar results: stars reddened above E(J — H) pa 
0.1 mag were well-fit by a reddening law of 
pa 1.20, while stars with lower reddening were 
consistent with noise. 

3.2. NICER Exinction Map 

Using the E(m\ 1 — m\ 2 )-Ay ratios derived 
above, an extinction map, shown in Figure [HI 
was constructed using the NICER method (see 
Lombardi & Alves 2001 for a detailed descrip- 
tion). The intrinsic colors and color disper- 
sions, used to calculate the color excess for in- 
dividual stars, were derived from the control 
fields (Equations |H and [5]). Visual extinctions 
were then calculated for each star, and the data 
were binned onto 4' pixel grids with 2' spac- 
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Fig. 10. — The pixel extinction distributions 
of the the sigma-clipped map (Figure [9]). In- 
dicating a population of zero-reddening stars 
is a Gaussian with o equal to the mean of the 
propagated measurment errors of the data. 

ing. At this resolution, we can resolve struc- 
ture down to ~ 58 pc. Pixels were required 
to contain at least three stars for subsequent 
calculations; those with N < 3 were given a 
"not- a- number" value and interpolated over in 
the extinction map. Out of a total of 44,590 
(245 x 182), there are only 6,089 such pixels, 
primarily in the periphery of the map where 
the number density is low. 

At our resolution, we found a mean extinc- 
tion (Ay) of 0.38 mag and standard deviation 
of 0.57 mag. The distribution of measured val- 
ues is displayed by the histogram in Figure [TOl 
"Negative" extinction pixels in Figure M corre- 
spond to regions that are bluer than the av- 
erage intrinsic colors or are due to errors in 
the measurements. For clarity, a map show- 
ing physical extinctions Ay > is shown in 
Figure [TTJ The most conspicuous feature of 
the map, which shows extinction due to both 
the atomic and molecular ISM of the LMC, is 

30 Doradus at a ~ 5 h 40 m and 5 69° to 

5 ~ -71° (J2000.0). Other filamentary and 
arc-like structures are seen, for instance, at 



a ~ 5 h 35 m and 5 ~ -68°20' to -69°20' The 
"CO arc" identified by Mizuno et al. fl200lf ) 



is also present in our map, extending south- 
east of 30 Doradus. Extending southwest from 
(a, 5) ~ (5 h ,-66°) to (4 h 40 m , -68°30') is a 
loose "string" of highly extincted clumps, some 
of which are associated with H I peaks. But 
there are no peaks in either the H I or CO maps 
corresponding to the rather large (~ 175 x 90 
pc) clump at (4 h 40 m , -68°30'). We do not find 
distinct, high- Ay features located toward the 
bar at (a, S) ~ (5 h 20 m , -70°). 

3.2.1. Uncertainties and noise 

The average extinction of a given pixel is 
estimated by the weighted mean of all stars 
falling within it: 



Ef = i WM 



(Ay) 



(7) 



where the weight for the ith star is given by 
Wi = 1/of, and of is the variance due to the 
uncertainties in photometry and intrinsic star 
colors, of = al hot t + o c 2 olor 4 . The error in the 
mean, or the observed variance, for the pixel is 
given by 



a' 



EllWi({Ay)-Ay^ 

(Ef=i wi) 



The signal-to-noise ratio (SNR) of the final 
map is dependent upon the smoothing tech- 
nique used. In order to eliminate suspected 
foreground stars that are expected to increase 
the uncertainty in Ay especially in highly ex- 
tincted regions, we would have liked to smooth 
the map using sigma-clipping. In this tech- 
nique, the local uncertainty in (Ay), o~a v , is 
estimated, and then Equation [7] is recalculated 
using only stars that fall within in-o^ of the 
first estimate of Ay. 

When the distribution of control field col- 
ors is sufficiently narrow, or when extinction 
caused by the clouds is sufficiently high and the 
fraction of foreground stars low, sigma-clipping 
is a practical solution to removing foreground 
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Fig. 11. — Sigma-clipped NICER extinction map of the LMC at a resolution of 4 arcmin. This 
map shows physically allowable extinctions, i.e., Ay > mag. 



stars that contribute to the noise but not the 
signal. Otherwise, especially when the latter 
case is untrue, there is the risk of underesti- 
mating the extinction. This is because sigma- 
clipping eliminates outliers with both high and 
low extinction. 

The trend of uncertainty increasing with ex- 
tinction (see Figure [T2|) results from a num- 
ber of sources of error, including the uncer- 
tainty of the intrinsic colors, small-scale struc- 
ture within pixels, photometric errors in the 



individual bands, and foreground stars within 
the LMC. The latter two uncertainties are ex- 
pected to increase in highly extincted regions. 
The stellar number density in high- Ay pixels is 
comparable to low- Ay areas of the map (~ 15 
stars per pixel and ~ 21 stars per pixel for 
Ay > 2 and Ay < 2 mag, respectively), and 
so Poisson noise is not the main contributor to 
the increased dispersion. 

In the NICER method, of is calculated for 
each star from the covariance matrices of the 
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Table 2: Summary of LMC reddenings. 
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Fig. 12. — The dispersion of measured extinc- 
tions function of Ay in bins of 0.10 mag. 
The black squares show cx mean , as calculated 
by Equation [7J averaged in each bin, with the 
relative scatters indicated by error bars. For 
comparison, we plot the mean error due only 
to scatter in intrinsic color and that due only 
to photometry, indicated by red and blue tri- 
angles, respectively. 

intrinsic scatter in colors and photometric er- 
rors. The corresponding measurement error of 
a pixel is given by a^ eas = 1/ W f- The 
following analysis of cr^ eas and the total disper- 
sion of a pixel, cr^ ean (Equation [S]), can help us 
determine which types of uncertainties domi- 
nate. 

As shown in Figure [12], in low-extinction 
regions of our map, measurement (i.e., pho- 
tometric and intrinsic color) errors contribute 
to a larger fraction of the overall uncertainty 
than in high-extinction regions. As the extinc- 
tion increases, the dispersion due to small-scale 
structure and foreground stars within the LMC 
increasingly dominate the uncertainty. As dis- 
cussed in previous studies (e.g. Lada et al. 
1994; Alves et al. 1998; Lombardi et al. 2006), 
these two phenomena are difficult to disen- 
tangle since they result in similar statistics. 
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a Except for the Isserstedt (1975) value, a maximum 
detection, values listed are means or ranges of E(B — V). 

b In converting our derived (Ay) to (E(B — V)), we used 
R = Ay / E(B — V) = 3.1, as did the authors listed here. 
The use of Gordon et al.'s more recent value calculated for 
the LMC, R = 3.41, would result in (E(B - V)) = 0.11 mag 

c These authors did not correct for Galactic foreground 
reddening. 

References. — (1) Feast et al. (1960); (2) Isserstedt 
(1975); (3) Isserstedt & Kohl (1984); (4) Grieve 
& Madore (1986); (5) Hill et al. (1994); (6) Massey 
et al. (1995); (7) Oestreicher & Schmidt-Kaler (1996); 

(8) Harris,Zaritsky & Thompson (1997) ; 

(9) Zaritsky (1999) 



For instance, an un-reddened star appearing 
through a "hole" in a GMC is indstinguish- 
able from a foreground star in front of the same 
cloud. 

We check the accuracy of our estimated un- 
certainties by examining the distribution of 
pixels with Ay < 0, The appearance of a signif- 
icant number (~ 1/5) of such pixels (see Fig- 
ure dCO) is consistent with a noise population 
with a standard deviation characterized by the 
mean measurement error. That is, stars with 
Ay < are well-fit by a Gaussian with a stan- 
dard deviation equal to the propagated mea- 
surement uncertainties, (<7 mea s) = 0.39 mag. 
Recall, the dispersion of the control group 
colors — the value that determines the signifi- 
cance level of the extinction measurements — 
is not much less than the dispersion of the 
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field star colors (see Figure [6] in § 12.21) . With 
deeper observations and better measurements 
of the intrinsic colors of a well-defined popula- 
tion, these statistics could be improved upon. 
Meanwhile, as indicated by Figure [TOj we are 
unable to distinguish real signal from noise be- 
low Ay « 1 mag. With an overall 1-a noise 
level of 0.57, the average SNR of the entire ex- 
tinction map is 0.66; it is 2.5 for those regions 
of the map where Ay > 1 mag. 

We note that the column densities derived in 
this paper may be systematically lowered if the 
extinction coefficients derived in § 13.11 are too 
large. This would be the case if we assumed too 
low of a gas-to-dust-ratio, (see Equation [6]). 

However, the measured columns are more 
likely to be lower limits, for the following rea- 
sons. (1) Foreground contamination results in 
the dilution of signal, and thus, lower values of 
Ay. (2) The most highly extincted stars likely 
went undetected. In § 3.3 we show that our 
measurements of Ay are relatively insensitive 
where Ico > 2 K km s _1 ; deeper photometry 
in these regions would raise the measured ex- 
tinction. 

3.2.2. Comparison with previous studies 

The mean internal extinction we derive here 
is in overall agreement with earlier measure- 
ments. Table [2] summarizes the mean redden- 
ings given by previous authors, each of whom 
used some form of the color excess technique — 
subtracting intrinsic from observed colors — to 
derive a reddening distribution. All used UBV 
photometry of early-type stars for their red- 
dening measurements. They were necessarily 
observing the most luminous stars, which may 
have additional circumstellar reddening due to 
heavy mass-loss. This could lead to overesti- 
mating the interstellar reddening, as Oestre- 
icher & Schmidt-Kaler (1996) point out. 

Other common sources of error and inaccu- 
racy in the results of previous authors included 

(i) the assumption of Galactic intrinsic colors, 

(ii) not accounting for the Galactic foreground 
contribution to reddening, (iii) not addressing 
selection effects related to completeness, (iv) 



small data sets, and (v) limited spatial cover- 
age. Adopting Galactic intrinsic colors (e.g., 
Feast et al. 1960, Isserstedt 1975, Isserstedt 
& Kohl 1984) leads to an underestimate of the 
reddening since MW supergiants (the stellar 
types used by these authors) are intriniscially 
redder in color than their LMC counterparts, 
especially in UBV photometry. Not account- 
ing for the completeness limit (e.g., Massey 
et al. 1995) also biases the measurements to- 
ward lower extinction, since highly reddened 
faint stars are underrepresented (see § 2.2.1). 
The extinction will be overestimated, however, 
if the Galactic foreground contribution to red- 
dening isn't corrected (e.g., Hill et al., Massey 
et al. 1995). 

Inaccuracies in the overall spatial distribu- 
tion of reddening may result from items (iv) 
and (v). Feast et al. (1960) had a small data 
set (108 supergiants) sampled from across the 
LMC. Hill et al. (1994) used 7 OB associa- 
tions, 3 of which are associated with 30 Do- 
radus, one of the most highly reddened regions 
of the galaxy. Harris et al. (1997) and Zaritsky 
(1999), as part of connected studies, produced 
maps of variable spatial resolution which cov- 
ered areas of ~ 2° x 2° and ~ 4° x 2.7°, respec- 
tively. The most comparable map to ours is 
that created by Oestreicher & Schmidt-Kaler 
(1996) from 1,507 O-A stars throughout the 
LMC. They did not grid their map, but rather 
analyzed the distribution of reddened stars. 
Like us, they found reddening throughout the 
LMC, including in peripheral areas not associ- 
ated with H I or CO peaks. They also did not 
find much reddening in the bar. 

3.3. The X-Factor 

The X-factor is defined as 

= ggg) [cm' 2 ] 
A " Icq [K km s- 1 ] ' [J) 

To determine A L MC) both the CO data 
and NICER extinction map were convolved 
to maps of the same pixel size (4') and grid 
spacing (2') and then compared on a pixel- 
by-pixel basis. These maps are displayed in 
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Fig. 13. — This map shows extinction due to H 2 . It is created by applying Equation [10] to the 
NICER extinction map (Figure [9]), using only pixels with tot al Ay > 1 mag. Overlaid are the 
CO contours of molecular clouds identified by Mizuno et al. (120011 ). 



Figure [131 For the reasons discussed in the 
previous section, only pixels in which Ay > 1 
mag (or, SNR> 2.5) are used. A linear rela- 
tionship between Ay and Iqo would allow us to 
derive X by fitting the data. However, above 
about 2 K km s" 1 , as shown in Figure [HI there 
does not appear to be any linear dependence 
between ico an d Ay. This lack of a linear cor- 
relation at high column densities most likely 
results from insufficient sensitivity in our ex- 



tinction measurements. The global mean SNR 
of the Ay map (0.66) surpasses that of the 
Ico map (0.23). But the comparative sen- 
sitivities of the two measurements appear to 
depend on the molecular content in a given 
area. In high column density regions, above 
Ico = 2 K km s _1 , these SNRs are 1.60 for 
Ay and 6.5 for Jco- I n low column density re- 
gions, Jco < 2 K km s _1 , the trend is reversed 
and the ratios are 0.82 and 0.18, respectively. 
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Fig. 14. — Plotted is the iccr^4v r elation for 
the analyzed regions (see the text for details). 
Data are binned along the Jco axis in bins of 
A/ co = 1.0 K km s _1 . 

In order to take full advantage of the data 
in the regions of most concern to us, i.e., re- 
gions of highest extinction, we proceed by de- 
riving the X-factor statistically. Of the 4,396 
pixels with Ay > 1 mag, there are 271 remain- 
ing corresponding to positions where Jco > 2 
K km s _1 . 

The total contribution of hydrogen gas to 
the column density, X(H tot ) = JV(H I) + 
2X(H 2 ), may be obtained from Ay. Statisti- 
cally, along any given line-of-sight toward the 
LMC, half the total column contributes to the 

measured extinction. That is, X(Ht ot v ' ) ) = 
iV(H tot )/2. In terms of the average gas-to -dust 
ratio calculated by Gordon et al. (120031 ). the 
column density of molecular hydrogen is 



X(H 2 )[cnT 



3.25 x 10 2 M v [mag 
iiV(Hl)[cm- 2 ]. 



(10) 



Figure [13] displays the map of visual extinc- 
tion due to H 2 , created using this equation. 
Finally, the mean of Equation M is calculated 
for pixels where Ay > 1 mag and ico — 2 



K km s _1 . We note that the suitability of 
this method relies on 2iV(H2) being signif- 
icantly larger than iV(H I) since, otherwise, 
large errors would result from the subtraction 
in Equation [TUJ For our measurements to- 
wards regions of the LMC where Ay > 1 mag, 
2X(H 2 )/X(H tot ) = 73%±17%, indicating that 
it is reasonable to proceed with Equation [101 

In units of 10 20 cm -2 (K km s -1 )" 1 , we 
found that Xlmc = 9.3 ± 0.4, where the er- 
ror around the mean (i.e., the random error) 
was adopted as the uncertainty. Values range 
from 0.10 to 47.1, and the l-er scatter is 6.9. 
Figure [TS] displays this distribution. In the 
same regions, we measure a total H 2 mass of 
(4.5 ±0.2) x 1O 7 M assuming a distance to the 
LMC of 50 kpc. 
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Fig. 15. — The distribution of Xlmc for the 
analyzed regions. The error bars are the devi- 
ation in the mean number of counts for each 
bin. 



3.3.1. Comparison with previous studies 

Our value of Xlmc, ~ 4.7Xmw, falls in 
between that of Mizuno et al. (2001) (~ 
4.5X M w), who calculated it from high resolu- 
tion ( 2/6) CO observations, and that of Israel 
( 119971 ) (~ 6.5X MW ; see Table [3]). The former 
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group assumed that CO clouds in the LMC 
are virialized and calculated the A-factor us- 
ing A = M vir /Lco, where M vir is the virial 
mass of a cloud and Leo is the CO luminos- 
ity. Application of the virial theorem to com- 
plex molecular structures, known to exist in 
the LMC, may be inappropriate. There is no 
a priori reason why the elongated structures 
associated with 30 Doradus and the supershell 
regions, for instance, should be virialized. Our 
value is also significantly larger than the Blitz 
et al. (2006), who did a revised calculation 
based on the Mizuno et al. (2001) observa- 
tions, of ~ 2.7A M \y- If the GMCs were not 
virialized, we would have expected in our cal- 
culations a smaller value of A L mc than that 
obtained by Blitz et al. (2006). Still, their def- 
inition of what constitutes a molecular cloud 
may have led to an underestimate in the total 
molecular mass. 

Table 3: Summary of LMC A-factors. 
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However, the molecular mass may be over- 
estimated when, at limited resolution, unre- 
lated clouds at different velocities are blended 
within a beam. At their resolution (8/8), Co- 



hen et al. (119881 ) were able to resolve only the 
largest cloud complexes (~ 140 pc). They de- 
rived a value of A L mc nearly 2 times our own, 
though their results did not assume virializa- 



tion. Finding that LMC clouds are roughly 
6 times fainter in C O than Galactic clouds, 
Cohen et al. (119881 ) reasoned that Almc ~ 
6Xmw , where they adopted Bloemen et al.'s 
(119861 ) value of A M \y = 2.8. But using the 
updated value of A MW = 2.0 (Bloemen 1989), 
would result in a value (A L mc = 12.1) a bit 
closer to ours. 

Garay et al. dl993h . using the same logic 
and value of A M \y as Cohen, found that 
A" T; m c ~ 3 9 in the region of 30 Doradus. Is- 
rael (119971 ) found an even g reater value (84), 
while Mizuno et al. (120011 ) observed no sig- 



nificant difference in the A-factor there com- 
pared to other regions. The highest value 
we measure in this region is 33.7, located 
at (a, 5) = (5 h 43 m 51 s ,-69°26 m ). Our rel- 
atively low values in this region may result 
from underestimating A(H 2 ) (see § 13. 3j) . The 
largest value we findoverall, 47.1, is located 
at (5 h 35 m 05 s , -67 d 38 m ) and is associated with 



a cloud identified by Mizuno et al. (120011 ) as 
having an integrated intensity of 4.2 K km s . 

Deriving the A-factor independently of viri- 
alization assumptions by calculating A(H2) 
from a comparison of f ar-IR surface bright- 
ness and A(H I), Israel (119971 ) also measured 
a larger value (~ 6.5A MW ) than our own. He 
gives a number of reasons for why his mea- 
surements of A(H 2 ), and therefore Almc, are 
lower limits. But it appears that his Ic q mea - 
surements, borrowed from Cohen et al. (119881 ). 
are the reason for his comparatively large value 
of Almc, since A oc I^X. The mean inte- 



grated intensity of the 22 clouds Israel ( 119971 ) 
uses to determine A L mc is 1-40 K km s _1 ; the 
mean I co detected in co rresponding positions 
by Mizuno et al. (l200lh is ~ 8.5 K km s _1 . 
This dimu nition in measured CO flux by Is- 
rael (119971 ) was likely a result of his convolving 
Cohen et al. 's (119881 ) original map to a coarser 
resolution of 15'. 

Previous authors' derivations of A L mc, 
whether over- or underestimated, necessarily 
tie to their measu rements of to tal H 2 mass. 
Cohen et al. (Il988|), Israel (Il997h . and Mizuno 
et al. (l200lh derive values of 1.4 x 10 8 M , 



16 



(1.0 ±0.3) x 10 s M , and 4-7xl0 7 M , respec- 
tively. Our value, (4.5 ±0.2) x 10 7 M , as well 
as our H2-to-H I mass ratios are consistent with 
Mizuno et al. f|200lh . Based on the H I data 
of Staveley-Smith et al. (2003), we measure a 
global ratio of ~ 10%. In high-column density 
regions, where iV(H I 



> 1.5 x 10 21 cm 2 as 



defined by Mizuno et al. (120011 ) . we find 
tio of ~ 50%. This consistency in mean global 
measurements is probably due to our relative 
sensitivities: the CO observations are more 
sensitive in high-column (i.e., high-mass) re- 
gions, whereas our NICER map is more sensi- 
tive in the more numerous, low-column regions 
of the LMC. 

Further in keeping with the relative sensi- 
tivities of the NICER technique versus CO ob- 
servations, discussed in § 3.3, is the appear- 
ance of our maps. With the exception of the 
30 Doradus region, there are few coincidences 
between the peaks of our H 2 extinctio n ma p 
and GMCs identified by Mizuno et al. tOOlh . 
as shown in Figure [13l The group catalogued 
a total of 107 molecular clouds, including 55 
"large" molecular clouds. But throughout the 
LMC, we measure many places where the H 2 is 
not associated with CO. In a few regions, ex- 
emplified by the areas (a, 5) ~ (5 h 40 m , -69°) 
and ~ (5 h 30 m ,-68°30'), H 2 and CO overlap, 
with the former extending further than the lat- 
ter. More commonly, CO and H 2 regions ap- 
pear adjacent to one another. This suggests 
that apparently discrete structures are some- 
times part of larger complexes that, observed 
with either method alone, cannot be seen in 
their entirety. For instance, the three clumps 
extending from ~ (5 h 25 m , —66°) — an H 2 peak 
bordered by two CO GMCs — may be part of 
one, elongated structure with varying column 
density. 

4. Conclusions 

We have used the NICE (Lada et al. 1994) 
and NICER (Lombardi & Alves 2001) tech- 
niques to derive, respectively, the NIR redden- 
ing law and Ay distribution of the LMC. De- 



spite the large number of assumptions and un- 
certainties necessarily introduced in applying 
these methods to extragalactic sources, results 
were largely consistent with previous studies. 
Moreover, with the unprecedentedly large data 
set from 2MASS, the measured Ay distribution 
may indeed be an improvement. Our main re- 
sults are as follows: 

1. By comparing the HI and E(m\ 1 — 
m\ 2 ) distributions (the latter derived 
with the NICE method) on a pixel- 
by-pixel basis, we calculated the NIR 
extinction coefficients for the LMC to 
be Ay/E(H - K) = 20.83 ± 0.52 and 
A Y /E(J - H) = 17.30 ± 0.46. Divid- 
ing these results E(J — H) / E(H — K) = 
1.20 ± 0.04 gives a result consistent with 
independently measured reddening laws. 

2. Using the NICER technique, we created 
an Ay map whose spatial distribution 
showed many areas of high column den- 
sity gas across the face of the LMC that 
were traced by little CO emission. The 
mean value of 0.38 mag is consistent with 
previous results. 

3. We created an v4 v (H 2 ) map and com- 
pared it to th e Icq map provided by 
Mizuno et al. (120011 ). In regions of sig- 
nificant CO emission we found that X- 
factorof (9.3±0.4)xl0 20 cm- 2 (K km s -1 ] 
and M(H 2 ) = (4.5 ± 0.2) x 10 7 M . ' 
Our value of Xlmc is less than some 
authors who, using other methods, prob- 
ably overestimated it. Since our value 
is greater than that calculated by both 
Mizun o et al. (2001) and Blitz et al. 
( 120061 ). this is consistent with their as- 



sumption the LMC GMCs are virialized. 
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